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Abstract: Hartree-Fock theory is supposed to yield a picture of atomic shells 



> : which may or may not be filled aceording to the atom's position in the periodic 

Q ' table. We prove that shells are always completely filled in an exact Hartree-Fock 

O \ calculation. Our theorem generalizes to any system having a two-body interaction 



that, like the Coulomb potential, is repulsive. 



^ ■ The Hartree-Fock (HF) variational calculation provides an approximate determination 

of ground states and ground state energies of quantum mechanical systems such as atoms 



^ \ and molecules, and is widely used in physics and chemistry. 



The picture of atoms it is supposed to yield is one of "shells" (or degenerate eigenvalues 
of the HF operator) which may or may not be filled according to the position of the atom 
in the periodic table. (Indeed, the concept of shells stems from HF theory itself, for shells 
are not overwhelmingly evident, or even precisely defined in the exact many-body wave 
function.) Despite years of attention to this subject, it is surprising that this notion of 
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unfilled shells has gone unchallenged. A very simple proof, which we give here, shows that 
shells are always filled in HF theory. In other words, the degeneracy of the highest filled 
level of the HF operator is always exactly what is needed to accommodate the assumed 
number of electrons. 

The theorem and its proof given below obviously generalize to any system in which 
the two-body interaction V is repulsive, i.e., positive definite as an operator on the two- 
particle Hilbert space. In particular, V is allowed to be spin dependent and to contain 
projection operators, as in the nuclear physics setting. The electronic Coulomb repulsion, 
for example, satisfies this positivity condition. The one-body part of the Hamiltonian can 
be arbitrary. For convenience and because of its familiarity, we use the atomic Hamiltonian 
as an illustration. 

Thus, we consider a Hamiltonian 



acting on A'"-electron wave functions, i.e. wave functions ^'(ri, ui; . . . ; rjv, crjv) that are 
antisymmetric with respect to interchanging {vi^ai) with (rj,(jj). In the example of an 
atom with nuclear charge Ze, U and V would be given by U{y) = —Ze^/r and V{r, r') = 



To obtain an approximate value for the ground state energy, Eq — min(\E'|i/|\E')/(\E'|\&), 
the HF calculation restricts attention to Slater determinants, i.e., wave functions of the 
form 



in which fi,...,fN are orthonormal functions of space and spin: {fi\fj) — Sij. The 
approximate ground state energy is then given by the HF-energy, which is defined to be 





$(ri, (7i; . . . ; riv, ajv) = (iV!)-'/'Det{/i(r,-, aj)} 



(1) 



Enp = min($|iy|$), $ has the form (1). 



(2) 
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Any minimizer, $hf, i-e., a determinantal function satisfying Euf = (^hfI-H^I^hf), is a 
HF ground state. It may not be unique. We remark that mathematical precision actually 
requires an "infimum" rather than "minimum" in (2) because a HF ground state may not 
exist. This will be the case, e.g., for an atom with N > 2Z + 1, i.e., a very negative ion 
[1]. For neutral or positively ionized atoms and molecules, however, it was proved in [2] 
that a HF ground state does exist and, at least in this case, the word "minimum" in (2) is 
justified. 

If a HF ground state does exist, it necessarily obeys the HF (or self-consistent field) 
eigenfunction equations 

h^Vk = SkVk (3) 

for all 1 < k < N, where /i$ is the one-body operator defined by its action on an arbitrary 
function of one space-spin variable by 

{hMr,<T)=l-^A + U{v) + I J2 El^.(r',r)|V(r,r')dV]/(r,<7) 

y T=±i j=i j 

N 

- E E^^(r'^) hPjir',r)fir',T)V{r,r')d'r', 

r=±lj = l 

and where '^i, ■ ■ ■ denote the special N orthonormal functions comprising the energy 
minimizing Slater determinant ^»hf- The eigenvalues, Sk, of give us some insight into 
the possible energy levels for binding an extra electron, but that is not our concern here. 

Theorem: Assume that V is positive definite, i.e., for every nonzero function i/j of 
two space-spin variables 



(4) 



Let <f be an eigenfunction of with eigenvalue s (i.e., h^<fi = e<f) that is orthogonal to 
the minimizing set ^i, . . . , Vat, i-e., {^\^k) = for all 1 < k < N . Then e > Sk for all 
l<k< N. 
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Before proving this theorem let us point out its main corollaries. First, it implies that 
the functions 'Pi, ■ ■ ■ ,'Pn comprising $hf occupy the N lowest energy levels of the 
reader may or may not find this surprising, but we point out that there is no proof of this 
assertion without the assumption that V > 0. Our main point, however, is the second 
implication which, indeed, is surprising: $hf does not leave any degenerate level unfilled. 
There is a gap because e > for all k = 1,2, . . . , N . 

Proof of the theorem: Assume £i < £2 < • • • < ^at- We shall derive a contradiction to 
the assumption that e < Sn- First, we introduce some more notation. Denote e by Sjv+i 
and <f by <fN+i- Further for all 1 < /c < + 1 and 1 < / < A + 1, we define 



a,a'=±l •' 

Notice that V^^fe = and V\.^i > if k ^ I since V is positive definite. 

Now let 6 be the Slater determinant built from <fii, . . . <Pn-i, Vn+i, as in (1). One 
easily checks that 




and 



N N 



($HF|iy|$HF) 



fc=l k,l=l 



N-1 N-1 N-1 



(mm 



,N+1 



k=l k,l=l 1=1 



and, for 1 < < A + 1, 



N 



hk + Yl = {'fik\huF\'fik) = £k- 



(5) 



1=1 



Notice that the term I = k m. the sum in (5) does not contribute since Vk^k = 0. 
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Since $hf is the HF ground state it follows that (^hfI-H^I^hf) < (^l-f^l^) and, thus, 

< - ($HF|^i"|$HF) 

N-l 

= hN+1 — hN + ^ {Vi^N+1 — Vi^n) 

1=1 

= Sat+i — £n — Vn,n+i < —Vn,n+i- 
The last inequality uses the assumption Sn+i < Sat, but we then have the contradiction 
< -Vn,n+i < 0. QED 

The proof does not give a rigorous estimate of the gap Sn+i — £n, but it does show 
that the gap is at least Vjv,Ar+i, which is usually not a tiny quantity. Thus, even an 
"approximate" degeneracy is unlikely. 
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